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ABSTRACT 


The nonlinear space-time neutron flux equation with 
negative prompt feedback and one-group delayed neutron is 
reduced by the use of a nonlinear transformation to a par- 
tial differential equation, in which the nonlinear term 
represents a small perturbation. The general procedure of 
solution for the resulting weakly nonlinear initial-boundary- 
value problem is then estahlished by means of the method of 
successive approximation. Convergence of the analytical 
solution is also discussed. The solutions to a slab reactor 
core and a cylindrical reactor core are investigated here. 
Asymptotic stable equilibrium states are derived from each 
of these solutions. The present results are consistent with 
those obtained from previous stability analysis for the 
generalized buckling K greater or less than uf, the lowest 


eigenvalue of the associated linear Helmholtz equation. 
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USE La TTS 


The stability of nonlinear reactor kinetics has been 
investigated by several authors. Kastenberg and Chambré [1] 
have established stability theorems for the nonlinear space- 
time reactor kinetics equations by the method of comparison 
functions. Garabedian and Lynch [2] and Scalettar [3] have 
treated the nonlinear equations by modal expansions in 
Space and then solved the resulting ordinary differential 
equations numerically. Kaplan, Margolis and Harris [4] have 
applied a pertubation technique as well as a modal expansion 
in space to obtain approximate solutions for nonseparable 
transients. Other numerical techniques were summarized by 
Kaplan, Henry, Margolis and Taylor [5]. 

In this thesis, the reduction of the nonlinear initial- 
boundary-value problem to a weakly nonlinear one is estab- 
lished by a nonlinear transformation. The transformed 
equation is now solved analytically by means of the method 
of successive approximations. As an application of this 
method, the asymptotic stable equilibrium states are 
discussed for an infinite slab reactor core and a finite 


cylindrical core. 





TIA ER UA LONE xd I EOUZAEIONS 


A. GENERAL EQUATIONS 

In this thesis, the reactor is assumed to be an one- 
velocity, bare, homogeneous and Newtonian cooling system 
with one group of delayed neutrons. The equations describ- 
ing the neutron and energy distributions and the concen- 
tration for one-group delayed neutron precursor are given 


by Meghreblian and Holmes [6]: 


1 (rt) = vV Pict) + z, [(- 9) k-]éc9 e ACC) (2.1) 
9C t) "SSMO pert) (22h 
ESO) = et, $(r.t) a ^ ut ,t) - T, (6, €) (213) 


where the resonance-escape and fast nonleakage probabilities 
to thermal are assumed to be unity. The houndary conditions 


of equations (2.1), (2.2) and (2.3) are given by: 


nt = Ce rE + t) o for tye (2.4) 


and the initial conditions are 


$ (ro) 


where L. indicates all points on the boundary of the reactor. 


£(£) ; C(r,o) 2 q(£) ; T(r,o) z fce) (ao) 


The nonlinear effects on the kinetics of a nuclear 
reactor are usually described by a single parameter oT the 


feedback coefficient of the reactor. The effects of 





temperature on the multiplication factor are more pronounced 
than the effects of temperature on the other reactor 
parameters. Small changes in temperature cause small changes 
in all the reactor parameters. However, the multiplication 
factor is very close to unity so that the quantity, 

(1-8)k-1 is very sensitive to even small changes ink. Hence, 
it is logical to assume that all the nonlinear coupling 
appears through this multiplication factor and the feedback 


equation is given by 
A(T) = ko + Xo MH (2.6) 


The dependence of the multiplication factor on the 
temperature and the appearance of neutron flux ó in equations 
(2.2) and (2.3) make the system defined by equations (2.1), 
(2.2) and (2.3) intrinsically nonlinear and coupled. 


Introducing the new dimensionless variables 





¿A t 
(2. 79 
- 
x — D r 
the above system of equations becomes: 
dix «Y g? ^ _pyk-1|b(%,7 PU X Zu ' 
2$(x, 0) « v$(c70 « [o- R66 0) « à. ceno) (2.8) 
ES -. Act BR h(x ,T) (2.9) 
Beet) = — ACA = P o 
x) _ ez xz) hà e (2.10) 
Aur) = ae ) PV Y [fes )- iC | 
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Defining the constants: 


_ Bk EC - 3 M. A Dan 
MESS E h= eer Co pe 


equations (2.8), (2.9) and (2.10) can be written as follows: 


Be. ) = VPET) V [o-m&-i] $6 à vacuo O) 
2C) EIA AE (2.13) 
Bac = V OK = 8 [TS To (x) ] (2.14) 
: DENS DE DUE 
From these equations, putting EE 0, ues 0 and ás 0 


vields the steady state equations: 


o = v*$,0) « [uo &,-1]$ 09 * va,Co (2%) (2.15) 
0 = -AC (4) + % (X) (2.16) 
O = 3, p x) cath 3, Io) To (2) ] (2.17) 


which are assumed to be true for T < 0. 
Subtracting equation (2.15) (2 Te) an ON 


equations (2.12), (2.13) and (2.14) respectively yields: 





LET) — V^y(Et). [0-1 Jy am 4 VA,5 (4,7) (2 EB 
Dar) < MRT L X 9 Cen) (2.19) 
EU) == QV) — %,6(X%: 7) (2.20) 


1 





where E | = (BR! 
and (Ex) = Qia) $2) 
51,1) = Clen) - G 0) (2.20) 


— 


Gar) = Tat) To) 
and k(T) is defined by the feedback equation (2.6) 
R(T) — Roe €(x.x) (2.22) 
The boundary conditions are given by: 
WE, 7) = 52,7) = O(4,,7) = 0 (2.230 


where X. designates the position vector on the boundary of 


the reactor and the initial conditions are: 
(4,0) = F(%) 
b(x,o)— o (2.24) 
O NO 


Equations (2.19) and (2.20) with initial conditions (2.24) 


can be solved to obtain: 


KLL) 


T 
SEAT) yf e p(s 110) dr! (2525) 
07) = Bf eB EMP yet )d a! (2.26) 
oO 


It is customary to rewrite equation (2.26) in another form: 


12 





c 
— 
Ux 
e 
rt 
| 


[ 946309 6,04 (2.27) 


| Ds (aa: ) 


where ER (tec) U e (2-28) 


The explicit form of the kernel (2.28) depends on the feed- 
back model under consideration [7]. Substituting equations 


8122 (2.25) and {2.27) intofequation (2.18) leads’ to: 
Mur) < VLS pe] e o parem (corde 


+ Rl dR ride (2.29) 


which describes the increase in neutron flux above its 
steady value for one-group delayed neutron and Neutonian 
cooling model. 

The boundary and initial conditions for equation (2.29) 


are given by equations (2.23) and (2.24) 


RD OE L O (2.30) 


y(x,0) = PCS (2.31) 


B. ASSUMPTIONS 

Before the analysis of the reactor may be undertaken, 
two further assumptions must be made concerning the feedback 
model. | 

The first assumption states that the reactor temperature 
rises instantaneously with the flux (i.e., the prompt feed- 
back model). The kernel K (x-«') in equation (2.28) becomes 


[7] 


LZ 





SIT SON 


where 6(T) is the Dirac delta function, hence, the temper- 
ature difference in equation (2.27) can be written as 


follows: 


~ 


T 


THI C 
o 


o 


or le ar (28535 


The second assumption states that the temperature 
coefficient of a reactor is negative. This is the case of 


negative feedback: 


2R (o 


= (2.34) 


Therefore, the feedback coefficient a in equation (2.22) 


is also negative: 


A CO - (2.35) 


Inserting equation (2.32) into equation (2.29) provides: 


dar) c V'yqos [0-1 vos * (1-8), v^) 


ev. e" aly ae’ (2.36) 

Let (1-8), - ! = K e 
(1-6) %0% = -x (x>0) (2.38) 

VAS, = B (2.39) 
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Equation (2.36) yields: 






(2.40) 


Y 








T ot) 
PET) = VPET AKPT) emo" ponds! 


Y 





and its initial-boundary conditions are given by equation 
(2.30) and (2.31). Equation (2.40) describes the increase 
of neutron flux above its steady value for a negative prompt 
feedback and one-group delayed neutron model, the solution 


O which will be established in the next section. 


15 





FI REDUCTION OF A NONLINEAR SYSTEM TOZAZLINFAR SYSTEM 


A. NONLINEAR TRANSFORMATION 


The Wilhelm's transformation [8] is considered: 


an ee P u 
YLE, LT e 2) | 
or U(x,T) = ee) (332) 


Bc ate a i yl "C ) 


EL - neutron flux difference Uí(x,tT) is non-negative, it 


is logical to assume that: 


W(x, 7) Z O fer x in D and T DYO 
hence, 


i- «(i-e )ucct) »9 (3/3) 


The nonlinear transformation (3.1) is generally valid, since 
provides a unique interrelation between (x,t) > 0 and 
u(x,t) 2 0; witehmus®iree irom singWlarities for any point 


Met). Differentiating equation (3.1) yields: 


KT 


EX fot pu ion} © (3.4) 
OT ES re 

KT Z KT, KT e 
V^ E e Vu a za(i-e Je (vu) (3&5) 


an-A] eze] 
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Substituting equations (3.1), (3.4) and (3.5) into equation 


(2.40) provides: 


KT 


T kt z 
Es + Ku) e -— e vu ne Zzx«(i- €. )€ (vu) 
2^ IIT I-20- e ul? efi- S > 4r 


KT hel 
Ul 


Ke wu 
"aen E 2 (1- e (e Jul" 


X. = Lr «T! i / 
«ele T (3.6) 
ia (1- © met (4,7) 
¿KT 
Dividing both sides by = P 5 and simplifying 
I - K Qef ui 
give: 
ae et Š un. 

A ES )u T[ aes jr! (3.7) 
pt (ree -e m RU e Per 
where Ke = K+ $0 (ION 


Let T) AE END VO) (3.9a) 
i 


~ e ex Ur m T! d i 
and D = B igual. U(x, 7’) AT (3.9b) 
2 A Z (Se THUS, T') 
Equation (3.7) becomes: 
OU (Xx, T) = SA Br D, (£,7) + yp NO (3.10) 


OT 
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The initial and boundary conditions are: 


ONS SP) D O us. TU 


MR B O) = F(x) (3.12) 


B. COMPARISON OF TERMS OF THE TRANSFORMED EQUATION 

Before solving equation (3.10) by successive approxima- 
tions method, it is necessary to impose some conditions on 
the terms ni and n» such that they are much smaller than 
the other terms in equation (3.10). A comparison of the 
terms of equation (3.6) indicates that the terns(2 e" and 
(v ^g" have the same order of magnitude as the term Kue^ ^ 
Therefore, the nonlinear terms n4 and n4 are negligibly 
Eun equation (3.10) for any time N < T< œ if and only 


N, 


if the ratios and 








are small compared to 











Kuef" Kuek™ 
ENS i.e., for x in D and 0 < T < œ: 





KT 2 
E Ixl e EU vu) 2 (3) 


ti __ 








«^e SR KS e")u |u 
and 
E (Tv) YI 
E |. " KT zn (3.14) 
Kue""| ^ EC we : «JJ -S(i- en eure) $ 








Since u and Vu are regular, as t tends to zero, the numer- 
ators of expressions (3.13) and (3.14) approach to zero, 


hence , 
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ire n, ————» O 
To | KU E 
B m N — O 


t. GC T 





For the time being, it is assumed that the conditions (3.13) 
and (3.14) are satisfied. They will be verified after the 


solution to equation (3.10) is obtained in the next section. 


C. LINEAR EQUATION TO BE SOLVED 


Let e Am (3.158 


Equation (3.10) becomes: 


E E Tut) + Gs ac m] (3.16) 
CUTE 6 (3-179 
ea) = F(x) (3.18) 


The transformed equation (3.16) now can be solved by the 
method of successive approximations. In this approach, the 
small nonlinear term G[t,u(x,T)] is ated as a perturba- 
tion. Hence, the initial-boundary-value problem is written 


as follows: 


25 s Vut Gi, uue] (3.19) 
U(x T) = 0 (3.20) 
I ES) (3.21) 
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where Os = O 


g^ T. 1) 
Wan uc zt yur) (ea, PR dr 030205 
C d Z Ki 0 ANN +8 "[ (1 > Jf Lo zu -* (i-e )u, (t (ET) 
thus, the reduction of the nonlinear initial-boundary-value 
problem, equation (2.40) to a linear one, equation (3.19) is 
established by the nonlinear transformation, provided that 


muceeconditions (3.13) and (3.14) are satisfied. 


20 





IV. ANALYSIS OF A SLAB REACTOR CORE 


A. SOLUTION 

Consider an infinite homogeneous slab [Fig. 1]. The 
medium is of infinite extent in both y and z directions and 
of width 2a in the x direction (a is the dimensionless 
extrapolated half-thickness of the slab); the origin of the 
X axis is placed at the center of the slab. These specifi- 
cations reduce the spatial dependence of the flux to 


variations in x alone. 


Figure 1. A Bare Slab Reactor. 


Therefore, equations (3.19), (3.20) and (3.21) become: 








LL: e ^t. 
22:6) — Lin + b 
U, (19,7) = O (4.1) 





oe »* 
=> æ 





where | 76 = 0 


and 


2 


I T N E "T 


avo 


© ST 


E 
M KT e Ce ae b= | eee ee 
Ta u] u 


As seen in Section III, equation (4.1) will be solved by the 


method of successive approximations where the zero? approx- 


imation corresponds C = 0. 
1. Zeroth Approximation ape = 0) 


The zeroth approximation uo Gr, x) is defined by the 


following equations: 





o = i> 

OT u 2x7 

rn We (4.3) 
U, (x, 0) = ex) 


Equation (4.3) is of the homogeneous, parabolic type that 
can be solved by either separation of variables or using 
the Green's function [9]. Throughout this work, the homo- 


geneous equations will be solved by separation of variables. 





aa Wir) = AYTO) sel 
2 1 
hence, a a S and Ge — 1p i 5 
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Substituting equations (4.5) into equation (4.3) and then 


dividing by XT yield: 
a m - < M (4.6) 


Pneresisıs the eigenvaluesto besdetermined from the homo- 
geneous boundary conditions: 

Zus a) L MG (4.7) 
Equations (4.6) can be written as 2 ordinary differential 


equations: 


X 2HBHgX = 0 (4.8) 


Le 


i a o (4.9) 
The solution to equation (4.8) is 


X (x) = C, cosux + C, sinux (4.10) 


where Ci and cj are 2 arbitrary constants. The boundary 


conditions (4.7) give 


C,cespua + GSinM&à = 0 (4.11) 
and C,cosHa ~ © SinMa = 0 (4.12) 


Subtracting equation (4.12) from (4.11) leads to 


ec, sinka = O 
This implies Co = 0 or sina = 0. 
xf sinta = 0, then La = nly n®=90,1,2,... and equation 
(4.11) provides 
C, Cosnit = O 
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or ce = O 
and the general solution (4.10) is of the form: 


CX) E C, Sin CT x Gt. 13) 


This solution is unacceptable since uc Gt, c) is non-negative 
lu uu -*x = ar hence, 


C, C MEO (4.14) 


and equation (4.11) gives 


Il 
O 


C, cos pa 


Since Ci # 0, then 


| 
o 


Cospa 


or the eigenvalues are defined by 


x = (n) Iu ORAE (4.15) 


When equations (4.14) and (4.15) are inserted into equation 
(4.10) the solution of equation (4.8) is obtained as follows: 
X (x) 2 C, cos CLEDRA" (4.16) 
Za 


The solution to equation (4.9) is 


272 
- (2n+') EC [v 
42 


MS EE CIA) 


where c, is an arbitrary constant. Combining equations 


3 
(4.16) and (4.17) gives the general solution to equation (4,3) 
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o HA x 
a T » BG cosp,r (4.18) 


nzo 


where „ are the eigenvalues given hy expression (4.15). 
The coefficients B. are determined from the initial condi- 


tion: uo (x,0) ZSF (x) 


u, (x, o) 


F(x) = » Bn cos un x (4.19) 


h=0 
The above equation indicates that the Bs should be the 
coefficients of the Fourier's series expansion of function 


Ao RAe., 


JA; if E Gespa dg (4.20) 


Replacing B. in equation (4.18) by its value in equation 


(4.20) yields 


ec E us a 
U (x,t) = » e rl F(5,) cos pen y, de | 


B ZQ a 
= 
or pL 
Since the series X e cosu, x cosy 8 is uniformly con- 
n=0 


Bement tor =ya)< x < a and 0 « x « efcf. appendix A] the 


interchange of signs È and J is possible, then, 


2 ec e 
U, (x. X) = / E2 tà e^ cast, x oput |i (4.21) 
-A nzo 
Let | | 
Qc ETE SERE ne | (ON 
/ j b», E COS MAX COS E 
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Equation (4.21) can be written in the form: 


(4.23) 





° 
Í G(x,5,x)F(5) dg 
"å 





where G(x,&,T) is usually called the Green's function [8] 
associated with equation (4.3) 
fi Fixstiappreximation 
The first approximation is the solution to the 


following equation: 


Qu, | 2'u = 
Sr u isn © G(x. 7) 
u (+ E 
(Fam) z (4.24) 
uU, (x,o) = Fc) b 
where] is obtained from equation (4.2) by making i = l 
and substituting equation (4.23) for uo (Qt, 1). The#sslutıon 


w 


of equation (4.24) is the superposition of two solutions, 


uo and v to the following equations: 


Que E chus (4.25) 








xU NES 2v = 2% 47) (4.26) 
U, (ta,7) = Q l v(ta,) = 0 
Uo(x,0) = F(x) v(x, o) = Q 
u, = UL + Y (4.27) 


Equation (4.25) is identical to equation (4.3), hence 


the solution uo (x, 1) is given by expressions (4.18) or (4.23). 
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The solution to equation (4.26) is now assumed to be: 


NAE c4 2 A y (7) 205 Mp X (4.28) 
ri =o 


the H's are determined by equation (4.15). The source term 


$4 G1) in equation (4.26) is then expanded in terms of the 


eigenfunctions, cosy, x: 


tu) > 2 É, (x) cos us x (4.29) 
t a 
where Br = Z) E, (4, x) esp, d (4.30) 
e 
En. 


Substitute equations (4.28) and (4.29) into equation (4.26) 


to obtain: 


oo zn oy 

e 2 1 ES 
> A (v)cospuX = - y un An TICOS X + >» 3 (T)EOS MX 
PLZ O n-o P Z O 


: œ n 
where AL ft) = : 
= . 
. 2 1 E 
or ys [Ån Ba Ani) ~ € (7) feos, x =o (4.31) 


nzo 
Since the eigenfunctions, cos y,x are linearly independent, 


equality (4.31) implies that: 


2 
Ant) + Ry An tt) a (4.32) 


| 
OS 
e 
c) 
A 


The initial condition is given by: 


257 





V(x, O) = 0 «= Y), A,(o) cosy, x 


which is equivalent to: 
Ay, Co) NO QAM zer (4.33) 
Un 
Multiplying equation (4.32) by e and integrating from 0 


to t yields the solution to this equation subject to the 


Militia! condition (4.33): 
ue K pat 
Am -e” |& (t)e dt (4.34) 
O 


where É 7 (t) is defined by expression (4.30) 


2. E a 
A 7 E -Pa Y 25 es E 
„Were [e Ja] % Gt) cost at 


ote E 
offe € (5 ,t)cosu,£ dE d (4.35) 


Or 


Inserting equation (4.35) into expression (4.28) leads to 


the solution of equation (4.26) 


o9 = a 2 
vem eS eme (fer hene 
Q -a 


nzo 
Since £1 Gt) is well-behaved rre EXEC) Ethesunirorm 
= E) 
convergence of the series PX e COSH X COSH É? 
n=0 


[cf. appendix A] allows the interchange of the symbols of 


summation and integration in the above equation: 
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T eg có ix 
vit) S J J E =) ic TLE ld (4. 
4 La nzo " 


From the definition of the Green's function, equation (4 


equation (4.36) can be written as 


a 
V (x, X) -| [ews rotg Ag dt (4 


o “oy 
As seen from equation (4.27), the first approximation is 


given by 


(4 





T d. 
y A woa f [965.7 6 (£, t)dýdt 
Q +a 


ar ‚eh Approximation 


36) 


‚228 


37) 


<38) 


It is noted that equation (4.1) for the ¡Eh approx- 


imation is similar to equation (4.24) for the first approx- 


imation with the same initial boundary conditions. 


Furthermore, the Green's function associated with a differ- 


ential equation only depends on the form of that equation 


and the boundary conditions [10]. Therefore, the solutions 


to all possible equations (4.1) with different ©, (x,t) can 


be deduced from expression (4.38) by replacing &, (x,t) for 


ti (x,T) 


(4 


LU 





where ~¢, (E,t) is given by equations (4.2). Once u, (x,T) 


EN 


evaluated from equation (4.39), the i approximation of 


29 


=» 


is 





the neutron flux difference will be calculated from the non- 


finegar transformation (3.1): 


I< E 
Thus, the analytical solution for an infinite slab reactor 
core with negative prompt feedback and one-group delayed 
neutron is found. 
At this point, it is necessary to verify the conditions 
(3.13) and (3.14) and the convergence of the sequence 
tu, (x,t) }. All these questions are discussed in the next 


subsections. 


Se VERIFICATION OF ASSUMPTIONS 

As shown in part III-B, the conditions (3.13) and (3.14) 
are justified at tT = 0 and in its neighborhood provided that 
the function u(x,t) and its derivatives == are continuous. 
These conditions will be verified for any T. 

For the sake of simplicity, the zeroth and the first 
approximations are used here to calculate the values of 


(x,t) and u(x,T) respectively. 


Expression (3.22) yields the value Or ce 


where 
KT oy 
7. NEM 2« (€ T (4.41) 
ki É(1- eu, ] 
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and 


t A Krr Mod eoo / 
| = Se )ve | f € Cio (x, v ) dt (4.42) 
= = b I- SE (i-e yu (xx) 


The fundamental mode of uo (x, T) is given by equation (4.18) 


H r 
AICA TI = pe costa, x 


and so 


a 
NR 
geom - Apel bons 
Let Aa, (x) = 


=~ 6, ald (4.43) 


The function a, (x) is positively small for -a < x < a and of 


order of = B. since à is very small (or order about mee 


) 


oco, x = 6, <x d 


(4.44) 
Expression (4.41) and (4.42) are written as: 
Z $e, 
Bez) -MMe*.. 
Im. Zag p ee Jsintpx (4.45) 
E p Kt ane TR a cog * P5 
wee abs eng CTE et aae 
t x 5 ro "TU, e a A ee 
= BBE j-a€ a2 [eos oX) L et ye’ 
EL E 6 [ peret e H a a o 
-u?1 
If the term a, e in equations 


(4.45) and (4.46) is neglec- 
ted for simplicity, it follows that: 


3ı 





(K- 2us)T 


BS re (4.47) 
4A N Kf + Ace En re | 
x KR Ratan) ai Me de’ 
NL | e e g, e I +acne T[ nn C (4.48) 
S o C ct 


m > 
a. Case K = LR 


Inequality (4.47) becomes: 

















a v 
a Het kad Tt 
i"' a Z4. me a Ie. Bo Me e. 
B KCi+ a, J = 
or 
ZRT 
uB a 2X Bo eC a 2« b. E AR 
kuet d E 
Assume Q = NL Bo = ee K = 0.09, 2e = 1079, then 
-2 E 
Tt. eso mee "EC 
KU GET 0.09 9 
or 
ii a | or an 6 
tt y Tt > 
mute st < 








In the same manner, expression (4.48) is equivalent to 


| T 
P MT Kg ART 7 
I s fee uE e pg Be To AR 
(14-04) : 
E -KoT ALT 
“8. e -y 2 QN i) < aa " 





3/2 





Replacing K by (K+ E in the above inequality yields: 


mix L£he” 

















and 
1. A (+ pa JT KT he let 
I2. Ó e = e 
Iq e^ S ZW AK 
or 
N e" 
2 es for any A >= 0 
ue” Ped 
x -8 -6 
ai o o , A = 10 and K = 0.09, then, 
N, 2 
IEEE TUS CA C. e 
| KU e < lo *x 0.09 9 í 
l 
un & | for any T, Of TS WA (4.50) 
ku ct 





Hence, inequalities (4.49) and (4.50) lead to the following 
conclusion: 
| % 


—" Oo I OST (4.51) 
KUO i 





This means that the nonlinear transformation (3.1) may be 


applicable in this case. 
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2 
< 
b. Case K Ho 


Inequality (4.47) can be rewritten as 


an 
2 2 (K-24, )T 


ULL o wee 














and so, 
| -p5 2 
je 2X ze a zap H, E 
alten] < ne imt ony 
Mere: = 0.09, K= 0.08, a = 10 --, B = 107°, A foe 
O O C 
then, 
i «2l to -~i 
| n. 4 AMOO O CI E ro 
uo (0.08)? 
or 
-— it ferany TZ o (4.52) 
rea oF* 
Similarly, inequality (4.48) turns out to be: 
ze cet c (16, - - HS 7 œ -KT R< HTS 
Ml < Be je ir’ cR eH usn 
BEE x - ESk: a ee 
T o O C O : ` 


it can be assumed that 
(k,- we ) < O 


then, inequality (4.53) yields: 


IR < Eee pe B EA me 
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Or 











2. 
—( 2K- a) Ge 
| Y 2 4 o EC. 
Ee x (ug) 
2 Ho 
PLI 2K - Ho 2.0, I dcl. "EDS 3^ ! it is obtained from the 


above inequality that: 


FOr Aall T > 0 


Ies | L k(us-K,) 


- 8 -4 
lo pa X16 


a LIS Ww: 
E | o. 08(0.01) 8 


Kuekt 





This implies that 





|- ne | LK | for all t > 0 (4.54) 
Kıte 


Expressions (4.52) and (4.54) prove that the value of 


| (x,t) | can be considered DE small compared to 
Kue 
U 
eT OC any 0 < T<ow if Cane K < u^. 


2 


c. Case K > U 
O 


It follows from inequality (4.47) that 


2 (ic- 2Mo I 





E < aasi u e for any T > 0 
hence, 

En, 4 2x Bout g Me™ ^ EP. for any © 20 
ku e*t K? K^ j S 
Let a = lon, B = Mv u^ = 0.09;7 K = 0.093 
Then, 
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i ~2] is en e 
| a | N GENOA VO x 0.00, 7 2.1 * 1o 


— 











kue“ CON 
l 
or es Z | for any ™ > 0 (4.55) 
200m o. 
To prove the ratio ENT is very small compared to one for 


any T > 0, it is divided in three different cases for the 


sake of simplicity: 





(e MIT : 
ie = Gf a, e << [ c.e, 
(K- ME )T 
1 le Na gee 
a, E e og. 
or z E nl K Ve — Tal 95A 
(K-p>) a Po 6.005 lo x1o'? 


c cC 760607 


Expression (4.48) yields the inequality: 


RAR L S 
LN < 6*5. e | E Mo) de! 
a 


or 

x -KoT Le SIE hs * -Hoz 
B, in | b, 
7 PI e (e e Eod 

o- Mo 2 2 

and 

No x — Kr MN 
| | a G ry x K CKo- u5) pon M 


Assume #@* = 10°, K = 0.093, K,- us ~ 0.003 
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Then, 








-8 P? 
iis A 10 ioe or any O (4.56) 
Ue 0. 095x0.00% 
(K- Mo YT : NE 
E > If a, E >> E 2 
T >> l ln K Ya "ctt 
K - MS Lo 


It is obtained from inequality (4.48) that: 


s a IA A bi CR: < Nec 
pm | x E Gace, <i e Ar’ 


t “Cc 2 
LRA TRR caps 2c F. 
e ade’ 


o 


a 
« Pe 


2 
Ln? (k-2Mo)T 
Bx O 
len: 














then, 
e Mr * 
| Tl s. 4 C Boae. L p bea for ang tz o 
EUo** arm Ee 
Dr 
m. -8 o 22] za 
er ¢ 10 “to x IO ~ 10 <& | for any ™ (4.57) 
RE (0-093) (15%) 
2 
(K- )'t 
5. HUE a, e Es El erde, 
~ l g kK EG 
t x a AA 
kK- M: => 
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chen, 


-KaT f (roo) a o 
n Cathe [nm ERa e 


2 


< af EE 
Ko- HS 


It is deduced from the above inequality that: 














i x 
2. 4 ( 2 um ag cor any T yo 
IST 2 \ x M5) 
kue wi o” 
1 -b -4 
or (UE e | |. (4.58) 
Iu e** 0.092 x0.003 





Therefore, inequalities (4.55), (4.56), (4.57) and (4.58) 
shows that the conditions (3.13) and (3.14) are satisfied 
any t > OF) mn Other words, for all cases’, it can be 


written as follows: 


E «c <ul for any T > 0 (4.59) 
Kue 


C. CONVERGENCE OF THE SOLUTION 


Consider the sequence of successive approximations 


fu, (x,t) |u, (x,T) is defined by equation (4.39)} (4.60) 
T a. 
Ea - | | G (x,£, *- S6, (s, €) dg dt (4.39) 
o -él 
ER CREE 


oG 
Where GE =) = r3 e SOS a 


Z O 
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T TIE L 2 1 EE 1s)continuous in -a < x, <a 
Ra TOL appendix A] and bounded. Thus, the con- 
vergence of the sequence tu, (x,t) } reduces to that of 


{¢. (x,t) } or" (je uns ec, 
1 
KT 


Ku.e 
1 


Since the ratio is very small compared to unity 


ber O < tT < we have: 


KT 


E Ne CL cc 


L 


EOr 0 s T « Tir the function u; (x,T) is analytically well- 
behaved, hence bounded by a number M and the inequality 


above can be written as 


KT 
<< KM 


E 


L 


which implies that ane 7908 should be bounded by a quantity 


ET 


In addition, equation (4.39) gives the following recurrence 


Formula: 


m 
UL. — t | [Sa 5. (5, 9- t; (5,€)|4r 4€ (4.61) 
a - å- 


Bso, Zee 


where theG.'s are evaluated from expression (4.2). Since 


-?l and 8* about 10^? 


aq is of order of magnitude about 10 
[c£. appendix B], the first term in the right hand side of 
equation (4.2) can be neglected compared to the second term, 


l.e., 
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2 T Er (Te ; ( / 
yl 323 e C LE (4.62) 
9 = i ae NES Xo 


L As eT a z li- a jus 


L 


It can be seen from equation (4.62) that: 
El » 0 


since the integrand under the integral is non-negative. 


The zeroth approximation yte SEGA E E 


Next, equation (4.61) gives: 


T nà 
U, -Uuo = | i G(x,&, t-€)Ó, (5, c) dE 4€ 
o - A, 


where SITAS O 


and b, eC. > O 


Therefore, we obtain 
or Ui, 2 ug (4.63) 


The values ORE corresponding to u and u, are given by 


equation (4.62): 


tC) 
x ppi (eT u jan (4.64) 
|- Z val Ex J T) 
S E » 
Cone cub) ccr. 
= le dI YE =—_—_—_ EH (4.65) 
^g Ki ( zn a Pu, Ge!) 


The values of uo and uy do not much change the integrals 


in equations (4.64) and (4.65) but do change appreciably 
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the factors before the integrals. Furthermore, uy > uo’ 


hence, 
¡A ce (4.66) 


Combining equation (4.61) and inequality (4.66) leads to 


the following inequality: 


t= a > Oo 


or "PD MCN (4.67) 


Following the same procedure, it can be shown that: 


LNE Asc E Re $761 (4.68) 
and u,« «€ uU € € ui $ Ua S (4.69) 


The sequence ma is bounded and monotonic increasing, 
so it converges [11]. Therefore the convergence of the 
sequence tu, } is proved. 

ASYMPTOTIC EQUILIBRIUM STATES FOR AN INFINITE SLAB 

REACTOR 

Stability is one of fundamental problems in nuclear 
reactor design. Some questions are raised here. Is the 
System stable if the steady-state flux $6 (X) is perturbed 
locally or if the generalized buckling K is changed by 
changes in the absorption cross section due to control-rod 
movement? So, this sub-section is devoted to the investi- 
gation of all asymptotic equilibrium states (if any) from 


the solution obtained in paragraph A of this section. 


4] 





l. %Zeroth Approximation 


The zeroth approximation of the increase of neutron 
flux above its steady value is given by equations (4.18) 


amc (4.40):; 


ec = 
-(MA-I)* 
De BAe i COS Mi X 


gcn = Sa =. ete tt 

LA An na | x 

IE 2 Ple ”- e GSM 
nzo 


(4.70) 


where the B,'S are defined by equation (4.20). 


0 La "2 
a. Case K = ue 7 (53) 
All (ue - K) are positive for n > Ó since 
Cm cL QC Hac MEA nn 


Therefore, all exponential terms except the term in (Y 7 K) 
will become negligibly small after sufficiently long time 


and equation (4.70) reduces to 


WY, (x,t) ————.  -Bef9sfa* |, as T (4.71) 


E z A Cos. x 
which describes the RH approximation for an asymptotic 
equilibrium state. 
For the case of no delayed neutron (8 = 0) and no feed- 


back (a = 0), expression (4.71) becomes: 


Pol ——— O E s (4.72) 
za 


This steady state is exactly the same as that obtained in 


Reactor Analysis [6]. 
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2 2 
b. Case Ho Sn EN « My 


All exponents (ue —KJewith 1 =eGedre positive, 
hence the corresponding terms in equation (4.100) decrease 
with time t except the term in (ue - K) which increases with 


time since ue - K is negative. Thus, as t goes to infinity, 


Uo 1) approaches the following limit: 
(Moe) T 
Y lx, Y) Bee C OSH 6x 
JS Bar (H E oe 
1 + E : d 
: x -CH= KIT 
Since ES c. C os aX >> | as T — ee 


the stable equilibrium is 


(Xe) ——> (4.73) 


x 


For the linear case, the reactor would be called supercriti- 
Gal since the flux ye xc) will tend to infinity as seen in 
equation (4.73) by putting 02 0. Again, this result is 
consistent with previous Reactor Analysis [6]. 


C. Case K< L 


All exponential terms in equation (4.70) decrease 
to zero as time T increases to infinity. Thus, the increase 
of neutron flux above its steady value is zero at the 


asymptotic equilibrium state for the nonlinear case. 
SS as T —- «€ (4.74) 


mires, Situation is called subcritical ina linear reactor. 
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AE ESUEADpOroxrimation 


The first approximation of the increase of neutron 
flux above its steady value is obtained by inserting 


equation (4.38) into equation (4.40) and putting i = l: 


KT KT Pa 
Y, c, r ) = e uox) TE "s Qo "9t Ge) dee 
) -5 (1 -e Ju te" ZU e Jos ,7-€)6 (&,Odx4e 


(4.75) 
where u_(x,T), G(x,E,T-t) and 2, (E,t) are defined by equa- 
tions (4.18), (4.22) and equation (4.2) respectively. On 


the other hand, G4 ($,t) could be written in the form: 


EA) = Ne) + N, et) (4.76) 


where n, (E,t) and al are determined from equations 
(4.41) and (4.42). 

Only the fundamental mode is concerned with in the 
following calculations. Therefore, the values of us Ert), 


Qu 
jpg and G(x,&,t-t) are given by: 


: TU 
WE) = Be 65 H, R (4.77) 
mites 
9 V. = - in 
q 0 - -8.HMe0 sinp.E (4.78) 
~My (T- +) 
G(x E T-E) = He COS COSMAS (4.79) 
KT T rà 
Let JE (1 — — | | C L EE, E de (4.80) 
e -à 
T rà 
Lon = f | Sætt (s EAC (4.81) 
S S 
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hence, the relationship between T, (x,T) and T,(x,T) is 


established, 


a) = e Wa: (4.82) 


and equation (4.75) becomes: 


K 


Zc 
PAT) = E Ct o( x, Bruder) (4.83) 
m E. in 
1-2 (1-2 JU T) ALC ET Go 


— 


Substituting equation (4.77) for uo (x,1) into equation (4.83) 


yields: 
Lhe 6) 7 
V oso) e co ts O (4.84) 


S 2 
P EL -(ga7k2t ex «TI 
E bs |e -Q |cosp,x 2 x» Eel 


Taking the limit of the above equation, as Tt tends to 


infinity, gives: 


Mo- ET 
EE UU Meme cospex tT. [85 
i -MIT MS- EOT 
TO T- \-=B.[e E Jeosuyx- EI, + SI 


It is apparent that the temporal behavior of V4 (x, 1) depends 


on the generalized buckling K. The limits of Ii and Tos 


C approaches infinity, are evaluated as follows: 


as 


T pad ra "D | 
lim f | Soret ire dede = lim \ | G(x, t, 1- €)£, (x, €) dedg 
o 


(> -4. TAM ES = 


(4.86) 


Thus, the problem is reduced to finding the value of the 


integral, 
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TE E for t large 


which is the sum of two integrals: 


q - T 
jJ 663,7) 9 eu S | G(x.&,T-t)m5(x,t)dt€ 
6 o 


in view of equation (4.76). 


Equations (4.41), (4.77), (4.78) and equation (4.79) 


give: 

-0 poe ie 
feu. &-c)n (5,24 « - 2a e Hee Eus Esin Mu (Let . 
Let A((x,€) 2 - == Bo He cos MoXCoSHoE sin wok (4.87) 
and ale) = E B¿cos Ho E (4.88) 


The above equation can be written in the form: 


HIE =(PE-K)€ 


T BT T zer etu o "e 
fen dt < ALS, Ele J dt (4.89) 
o o mace? SACIE 


Similarly, taking account of equations (4.42), (4.77), (4.78) 


and equation (4.79) yields: 


"Qi-Ot Lute Holen? Sc nd D. de 
Jon, dt = = GNO p n je" [i- aye +a,€ j a SERV 
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Changing the order of integration in the above equation and 


letting 
Ve A 
Pu (x E) = £ E, cos, x cos" HE (4.90) 


leads to the following equation: 


T ST + 3 Ir POTE ye te c 
o o o e E K E 2 / 
en'dt- A,e H |. a ae ee 
J Te S TER IT) | Lı-a,e +2, ] (4.91) 


a. Case K =e 


As t is sufficiently large, the integrand of the 
integral in the right hand side of equation (4.89) is equiv- 


alent to: 


Br re aoe 
Se A.E = (4.92) 


— Em \ 
lz ae Mo act: \+ a (5) 


Substituting equation (4.92) into equation (4.89) gives: 


as. 


T - T 
fen diu Lc Atel? dt 
ô — o 14 45) 
pr. 
T j SM o TE 
Jena = A^ 3 (x, ES e T (4.93) 
where T | (4.94) 
| + e, (X) 


In the same manner, equation (4.91) can be written as 


follows: 
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C za mp e eor" "€ € 
-HoT Mo 2 (M,- X. i 
fenat K 16 | laes e Atdr 
ô 


o 1+4,(€) a 


when t is sufficiently large or after integrating: 


T (mc) EHE 
l enl d ~ air - — i mue jo (4.95) 
o Ko - “y Ko Sa H So 


But equation (3.8) is recalled: 


Ko = K+ AQ 
or oe a Ne 
since K- u = O 
hence, peri m X. (4.96) 


Inserting equation (4.96) into equation (4.95) and putting 
Aa ©) ae Fle) (4.97) 
4 A 


c 


leads to the following result: 
: 2 
T — Ho T Ae 
e ave 
Jj en ae SM Aa C (tT + = KL 
O 
Or, as T is sufficiently large: 
EE 


Tt 
vo 


Now, adding equation (4.93) to equation (4.98) yields: 
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r 
^ 


T T "C E 
fed at = Jenide+ [onde ES (A, +Aq je T (4.99) 
e o ò 


Integrating equation (4.99) with respect to & from -a to a 


will give the integral I,(x,t) as defined by equation (4.81) 


a ^t E 
L, (4 0) = fos, dede 2 Aste cx (4.100) 
-a o 
a 
where G (x) = IIT Ag (4.101) 
— d} 


The integral I, (x,t) is evaluated by using relation (4.82) 
and equation (4.100) 
KT - Hot 


kT 
Mer) eT ear) ae e Aw) T 


Since K = u then 
IS © (4.102) 


Replace equations (4.100) and (4.102) for I (X.T) and 


Ij (x,1) in equation (4.85) with K - ue to obtain: 


E. P Cx) = lim - Oe COS HOR SEY Se) T 

a ov Taan L+ ax) ana De aay Arg T x + SALT 
Since the term A, (x)T is much larger than the other terms 
in the right hand side of the above equation for tw large 
and A, (x) is always positive [cf. appendix B] for -a < x < a, 


Ehe limit of Y (X.T) becomes: 


Lim V, Cx, 1) = Lim pz 


Te ap Miño” 


49 





or Lim O a (4.103) 


«€ —» oo 


AMI V4 (x,) - = is the first approximation for a stable 
equilibrium state when the generalized buckling K is equal 


to the lowest eigenvalue, ue, of equation (4.3). 
2 
b. Case K^» U 
O 


If (K - uf) is positive, so also (Ko - u^) is, 
since K =K+A. 
O G 


For t large, the integrand under integral (4.89) 


becomes: 
EM) E nE PU one 
pa’. ea 3 gt = u (4.104) 
ut MI) E A X 
league Face I<) ACE) e es Kt d. (E) 
Similarly, for the integrands in equation (4.91) 
-(Ho-ko)t' (ude)! Ä 
e. d _ e Fonie a) 
- ur! zt EE — — 5. —— 
l-a x)£ La cg g nee AE) 
aad 
A (4.105) 
4 CE) 
and 
L 2 2 E -2(po-Kk)t 
(He- Kole te uote CHo 2 
P ahe cae | ze ate 
2 
2 (len Mer Ac VE 
Ela Ue i (4.106) 


Substituting equation (4.104) back into equation (4.89) 


gives: 
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z T 
T Hot NU 
fenas =x Aj £)€ ect m AN am (4.107) 
5 Q0 0 - 
where Acha) = Ale (4.108) 


e 0) 


Inserting equations (4.105) and (4.106) into equation (4.91) 


yields the following equation: 


E Buona dou e (He po ADE 
Gn dt Ae e a (ye Atd 
MEL E Lr 





rei, L 


then, after integrating and some manipulations, it is 


obtained that: 


^t EG z -(&- "x 
(«- 2k, )* AT Ho 
jet x Aaa 2 A mde tne | 
Ô 


(k-pž- 1.) ME x K- RE (4.109) 


The exponential terms in the bracket will approach zero as 


t tends to infinity. Therefore, equation (4.109) turns out 


to be: 
Y r 
fende s Be e $ = | 
E. K-po- A, Ae [Ss 
i T Br zu EE 
or fer, dt x P^ 2 
o + Ne (Ko) 
Let A;(x.€) = A1018) ute) (4.110) 
Ae (k= He ) 


It follows that: 


E (k - t4» )5 
Jon, de ~ Aj (x.5)e pa 
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T 
Again, the value of integral f G2,dt is obtained by summing 


O 
up equation (4.107) and equation (4.111): 


n 


* is © (- Ms D^ 
f Gt, dt Kc (x,£)€ 


¿2 t+ Ay | (4.112) 


When t is sufficiently large, the first term in the bracket 
is negligibly small, compared to the second term since 
ET — u^) is positive, hence 


(K-2ph )T 


{etd = Ade (4.113) 


The value of I, (x,T) is deduced from integrating equation 


64,113) with respect to &E from -a to a: 


E GR au DE 
mM.) 3 f [ot deaz = A 0e (4.114) 
-A o 
: a 
Where E coo J Ao (x, &)AE 


CT TES 
- A 


A. (x,T) is defined by equations (4.110), (4.90) and equa- 
tion (4.88). 


The value of I (X.T) is determined from relation (4.82) 


| rK HA 1 
S ~ AÅ (xn E E 


5 (4.116) 


Inserting the values of I,(x,T) and Ii (x,1) from equations 


(4.114) and (4.116) into equation (4.85) yields: 
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K- HEIT 2(X-m¿ Dx 


Lim W Cx, x) = Lim Boe eosu,x tAgixre 
= -p K- Ia Z TOP -ut 
oe scr Veneer ena eo a c MA, ez ctae 
or 


-(K-Ne)T 


1 
20«- p, Y 
Mos P (xv) Cis e [6,2 Coz Hah Ach 


= tpi) e -20k R3) LLK MZ) Y (pie -KT 
C2 o y [e i -AE p +A eL ~ Age +% Ag] 


T- 09 
K 


Te eT 
L me | Boe COSM, X + Ag 
— - ETERZE! —(2tc-MEDT (KE YT = 
LS LO H -42 +4, - E Age + =Ag 


Since (K - We) is positive, all exponential terms in 
the bracket go to zero as t approaches infinity. On the 
other hand, Ag (x) is always positive [cf. appendix B] for 


-a < x < a, thus, the limit of 94 (x,1) ter K > u^ "WIE 


MS UU OUR — IKE (4 IE 
T-> 09 = 
which describes the asymptotic equilibrium state of the 
increase of neutron flux, $(x,T) - $o (X). This is exactly 
the same result as obtained from equation (4.73) for the 


zeroth approximation. - 
2 
c. Case K < Ho 


The same procedure as carried out in parts (a) 
and (b), the stable equilibrium state, will be found by 


using equations (4.85), (4.89) and equation (4.91). 


L 2 
—( Hoy -KkYÉ — Mo E -kt 
n a e. (4.118) 


| ~= 


= c -(ut-x 2_ -KE 
pes A, e Ps = ae HM, )t go KE A e E A, 
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The above expression is obtained by dividing the numerator 


- (M5-K)t 
and the denominator of the left hand side by e : 


Since u^ - K > 0, equation (4.118) becomes for t large: 


Ue -me _(-k)e 
¿A e | Er ) (4.119) 


en 
m 

met _ zim-e)e Cos 28 
kace * ae € 


— 
— 


hence, integral (4.119) can he evaluated as follows: 


x TE -MoT a 
fon! de = e je Lee Ame [i-e | 
e © Bt oo e 
ct E 
Or fons ENE (4,120) 
Q 
where Ag (5) = me] (A, 12%) 


o-K 


Furthermore, integral (4.91) is equivalent to 


= C 
T -— U, T == es r TE / 
T eu) d ZEA e Ñ [e^ ef ett et dx 
6 : o mt 
2 L = 
-uT (Ma - Ko 
SER I i _ | 
U ere Mi Bas o 


It is assumed that: 


2. 


A 
then, for t large: 
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= 


wen 


re 
fen. K (ut = Kt 


R 


Aro RE E (4.122) 


n 


where A Ons) = Aas 8) (4.123) 
(Ho -Ke)” 


Adding equation (4.120) to equation (4.122) yields: 


x 2 
y -Ko T IIS 
GE, de = AQ0.*)Q0 A0, y)e 
Since ra > Keo 
thus, 
y — Kot 
J st, de Z Ajo (% B)e (4.124) 
and 
a T Lk. 
- Ex, ) zu [e£ 24s ze C. (4.125) 
= 6 
a 
where An) = | AE) AR (4,126) 
-à 
and 
- Ko) T = Act 
Ir) = re = T E Z T Cx) € i T 


Finally, inserting the values of I, (x,T) and I, (x,T) from 


2 
equation (4.125) and (4.127) into equation (4.85) yields the 


First approximation of the asymptotic equilibrium state: 
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Li L —(Ho-k)T TAE 
s zy sd - - UST LH E KT Kot IT 
To yo ! wa ns ^ axe = AWE + ZA 002 c 


Since all exponential terms decrease as T increases, then: 


Ere (x,t) = O (4.128) 


T> eo 


Thus, the same asymptotic state as for the zeroth approxima- 
tion Is obtained. The same conclusion can be made even if 
2 


Ko is assumed to be greater than ue, instead of Ko < E 


C PHYSICAL INTERPRETATION OF RESULTS 

It can be seen from the zeroth and first approximation 
that higher approximations are not required for asymptotic 
equilibrium states in a nonlinear system with negative 
prompt feedhack and one-group delayed neutron. The stable 
equilibrium states in the first approximation are verified 
by the results of Kastenberg and Chambre's stability 
analysis [1]. Furthermore, in the absence of feedback, the 
results obtained here reduce to those of the linear reactor 
analysis as mI EE part D-l. 

The equilibrium states are established in a nonlinear 
nuclear reactor at different levels which depend upon the 
value of the generalized buckling K. The results are 
interpreted as follows. 


For t < 0, the reactor was in steady-state with flux 


$o (X) : 


56 





At t = 0, the flux is raised to an amount of F(x) from 
an external source. What happens to the flux distribution 
after some time (t » 0)? 

IL NUS a the neutron flux will zeachkthe ¡stable 
equilibrium state of value (x) = $ 6 (X) + = after the 
transient given by equations (4.39) and (4.40) has subsided. 

g. FIER < y, thei neutron flux will return to the 
initial steady value $c (x) as equilibrium state. 

On the other hand, the equilibrium state is weakly 
dependent of delayed neutrons since the fraction of delayed 
neutrons (8) is very small compared to unity in the expres- 


sion of parameter a: 


c = 0 goo Y; | 
If delayed neutrons are neglected (8* = 0), equation 
(2.40) becomes 
p 2 a 
QU xs) = DH x ey Hie PEK ec Wa (4.129) 
oT x 


Hor K > (323 5, theorme I of Kastenberg and Chambré [1] gives 
the asymptotic stable equilibrium state of the system 


(4.129), which is the positive (nontrivial) solution of 


equation: 
lim [ 27 a e po) a, m] — 0 (4,130) 
> O ax? 

It is apparent that w(x) = = satisfies this equation, thus 

is a positive solution to equation (4.130). Therefore, the 
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result obtained in part 1) is verified by this theorem for 
K > (57 A Similarly, theorem II of Kastenberg and Chambré 
shows that the equilibrium state for K < (2^ is w(x) = 0; 
i.e., the neutron ES returns to the initial steady state 
$c (x). This is again consistent with the result in part 2). 
Theorem III states that if K - te, the asymptotical 
stable equilibrium of the system ahove is also the zero 
state (W(x) = 0). This completely differs from the present 
result as described in part 1). From the physical point of 
view, if the reactor eigenvalue is not changed, then the 


perturbed flux should eventually return to the initial 


equilibrium state W(x) = 0. The fact that for K = "d 


o* 


the first approximation contains a secular term at large 
t [cf. equation (4.102] indicates that this is a special 
case for which the Wilhelm's transformation may not be valid. 


To draw a definite conclusion at this point, it appears that 


2 


the first approximation of the solution for the case K = Ho 


requires further considerations. 
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Wome ANADYSTOSZOF A FINILIEJCYLINDKICAL REACTOR CORE 


A. SOLUTION 

Another geometric configuration which is of practical 
importance in Reactor Design is the cylinder of finite 
length. Consider therefore a bare cylinder of dimensionless 
extrapolated radius R and height 2%, as shown in Fig. 5.1, 
with the z direction along the axis of the cylinder and the 
origin of the coordinate system located at the center. The 
flux in this reactor depends upon r,z (cylindrical coordi- 
nates) and tT. The reactor equation is now [cf. equation 


(2.40) ] 


DY Z ET (5.1) 
EIU 7) = V o + K qp - cx wv - & fe pr, è, x’) dr’ $ 
C 
o 

Tor 62 ro cig 

where 
2 
mor < ro = (5.2) 





Figure 5.1. A Bare Cylindrical Reactor. 
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The initial and boundary conditions are: 


WCC) ee oy = FCT es) 


P(R,2,T) 


(IS) 
Y(0,2,T1) (5 Finite 


pcr, + 2,7) = O 


Using the same transformation (3.1), equations (5.1) and 
(5.3) are transformed to linear equations as follows 


E xeguation (3.31)1: 


PEL Ce! c Y S Via te Er) + t. On (5.4) 
OT 
Rr, o) Em 2 
DIR + me 
(555) 


U;(0,2,7) es finite 
Uml, t)= 0 


where A, = Q (5.6) 


and 


N if T -Koct-1!] 
2. za GR GF NES Zu] 2% usina rd (5.7) 
M = 5 t-l 
o 


K(i- Éc-e yu, , ] Aue U, K 
The methods used in the analysis of the infinite slab are 
easily extended to the analysis of the finite-cylinder 
reactor. The procedure is straight forward and involves the 


use of the technique of successive approximations and separa- 


ion of variables: 
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1. Zeroth Approximation te, 3 


The appropriate form of the transformed equation 


which gives the zeroth approximation is: 











2 
= al (ru) > Ze (5.8) 
eat r or oan o. 


= 
o 

co 
“4 
e» 
o 
Nat 

| 


Fecr, e) 


ECO Memo) 


Ô 


uUudlo, 2. T) CS finite (5.9) 


or, AS) 


l 
o 


A solution of this equation is assumed to be in the form: 


mo (tt, t) = Xn) Zar,» (5.10) 


The substitution of this expression into equation (5.8) 








yields: 
Tin t d ; 2” : (5.19 
— — Be - SNR - 
E C) CXe) ak X) BE PEN n 


where V is a constant to be determined from the boundary 


condition (5.9) or 


NN (5.12) 


It follows from equation (5.11) that: 


S" +U S 0 | (5.13) 
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and 


IAEA II y 
TÉ zx AEI ur. (5.14) 


y is another eigenvalue to be determined from the boundary 


conditions. 


SS) = Q ae 


X (0) is Finite 


Equation (5.14) can be written in 2 different equations: 


de(rEx') rE 20 (5.16) 
and 
Dan ja = 0 (5.17) 
Solving equations (5.12) and (5.13) gives: 
Hn = (2n +15 SS) 
and Zz) = C cosyn? (5.19) 


Equation (5.16) is a Bessel's equation, the solution of 


nich is of the form: 


X(r) = CI,(tr) + DY, (rr) (5.20) 


where C,D are 2 constants 


Xi) 2 on T 
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hence, D = 0 Since yg (0) — o and X(0) should be finite 
E ousudtryonct5sl5)].— Furthermore? the boundary condition 


X(R) = 0, yields: 
J, (7 By) - eG (5.21) 


If Qo, Yirg»: qr are the positive roots of equation 
(5.21), the solutions to equation (5.16) are given by equa- 


bison (5.20) 
rer ec, mr) (5.22) 


and the solutions to equation (5.17) are 
en EOS 
Tor) = E e (5.23) 
where E is any constant. 
Combining the product of solutions (5.19), (5.22) and 
(5,23) leads to the solution of equation (5.8) 
apo JT 


y 
U (rat) => A, e J (Y, r)cosh, + (5.24) 


m,nzo 


The constants Dee are defined by the initial condition 
r 
co of} 
EB, 2,0) = F(r,%) = »» Amn Jo mr) |cospn 2 


PL Z O rn 20 


This equation implies that: 


e3 2 
EJ. ar) = y [Ecos rcosunz ds = I D (5.2.5) 
ri Z G -2 
Mm = 0,1,2,--- 


) 
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Expression (5.25) is known as the Fourier-Bessel series 
representation of F (r). S are the positive roots of 


equation (5.21), the coefficients An , should have the 
r 


form [12] 
ray = — yee a. (an) 
a DEAN 
Or 
Z S : 
Eo V cee 13,0 n)FQ, Skssu,edpan (5.26) 
ART AR) \ J: (0,1) 11%) H t E A 


Substitution of expression (5.26) into equation (5.24) gives 


the solution uo (2,1): 
f Ye 
U Cn e, v) -[ | FAE) E de CHE Der) 1. Meeepn cp s Hin 5 27) 
o “E ETE J: (Win A) 


Let 


EY alts 
G(r,2,7,5,7t) = an Dé N cH G 1) cos, Beosp,, & (5.28) 


G(r,2,n,&,T) is called the Green's function of equation (5.8). 


Then, equation (5.27) becomes 


G(r, 2,» & ,v)F(, )ds an c 





2. First Approximation 


The first approximation satisfies the following 


equation: 
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ou $ 
pell dd eR 2 + Eres) E 
UL CC, 2, o) CoP) 
U,(R,2,7) E O 


(5.351 


O art ae te. 
U, (nd, 7) = O 


where €, (r,z,1) is determined from equation (5.7) hy making 
51. 
The solution of equation (5.20) is the sum of the solu- 


Exon u to equation (5.8) and the solution v to equations: 


V > OV D’ V 
VY (G2 05 = O 
AE, ) = O 
Ses 
v(o,2,7) is inite y 
VA 0 
u = Un + V | (5.34) 


The general solution to equatoin (5.32) is of the form: 


OS e (RES) 


zo 


2,3) = 


¿IN 
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Expanding t, muc bcrmswot Jo Gou. 2 yields: 


EE CY O roce, (5.36) 


nzo o 


oo 
where P» C J ("7 r) are the Fourier"'s coefficients of the 
mg mno m 


above series, i.e., 
o 


L 
D c A O nest 4s 


mr 6 
-£ 


This equation implies that C n n (T) are the Bessel's series 
, 


coefficients: 


ge 4 
Eo: zar YE (n, &,v)eesp,t AE A (5.37) 
0 J TAANE U SS 151 


Taking derivatives of equation (5.35) gives: 





d Bman Ji. GR c) cos M. Z- 


8 
I 
DI 


AT 
PR (TS D 
Ov = M n 
ar m m Bm, n J, Or) es um 
manso (5.38) 
92v = 2 
2p? ps 77 p E. PNE TS FA Mr) | cos pa, + 
rinh0 
9?v = 1 
ga ~ DI Hn Brin Jo (Imr) Sp e 
mm, ns 0 


Insert equations (5.38) and (5.36) into equation (5.32) to 


obtain: 
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m.nszo 


od 
DE LE DNA P, Cen | cosi. s 


ERr-Ecquation js trus for anv 0s xr siRiand. -L T < L if 


Endgeor TE: 





ab 2 L 
mr + ee = C ae (5.39) 


The initial conditions are given by: 


v(m #:0) = O = > Bm, (8) J, (Imr ) cosu, © 
or [o co) = 0 (5.40) 


The solution to equation (5.39) can be written as: 


= T e 
c FT (Tm + Mn Je 


© 


where bD. E are constants determined from the initial condi- 
, 


tions (5.40) 


B, M eo) = Do = WO 


i i 


Therefore, it is deduced from the above equations that 


T 
— vo E n S (T- R 
( e T 


E fe A E 


ð 


Replacing equation (5.37) for Cn n (E) leads to 
f 
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IL 


D R_ l J? 
-(U up 
2 ga C f. Ce C Rb o c AT CIA 


Then, it follows from equations (5.35) and (5.41) that 
(TE Dre 
V (S.Y) a f f Cn, pc) m 7; ye L 3 J. (Tm D (V. n)cosu »Ecosu, g Ag dnde 
min =o J,” (FR) 


or, in terms of the Green's function (5.28), the above 


expression becomes: 
rer pl 

Y(r,2,7) | \ [EBE Eleng mt) dg d de (5.42) 
6 e 7-4 


It is noted that this form of solution is analogous to that 


one derived in equation (4.37) for a slab reactor. 


Now, combining equations (5.29) and (5.42) yields the 


first approximation: 


(5.43) 





tT rR pt 
Ui (r,a,7t) U(r, xT) +f | [66e Eln., ded de 
. o o “—£ 


3. ith Approximation 
Since equation (5.4) is similar to equation (5.30) 
with the same initial boundary conditions, the ith approx- 
imation which is the solution to equation (5.4), can be 
obtained from expression (5.43) by substituting <, (n,&,t) 
for 27 (Mr8,t): 
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+4 
[ewan E, 7-)8, (1,5 6) dz dn de (5.44) 
J -4 





T PE 
tt: (C, 2,7) Bee, ker) «| | 
0 —o 


Where ~, (n,&,t) is defined by equation (5.7). Equations 
(5.44) and (5.7) will give all possible successive approxima- 
tions of the solution to equation (5.4) provided that 4, = 0. 

Finally, the ith approximation of the increase of neutron 
flux above its steady value is calculated from the transfor- 


mation (3.1) 


KT 
p rue) e m eura) (5.45) 
ex KT 
I\- £(1-e eer eco» 

Thus, the analytical solution for a finite-cylinder reactor 
core with negative prompt feedback and one-group delayed 
neutron is established. The convergence of the series 
tu; (r,z,1)] can be shown by the same method as indicated 


in part IV-C for a slab reactor. i 


B. DISCUSSION OF RESULTS 
From the on (5.44) and the transformation (5.45), 
the asymptotic stable equilibrium states can be found at 
different levels by the same approach as applied in a slab 
reactor. The proof is not repeated here but only results 
are summarized as follows. 
a. gf K = ue +e 
After a long time, the increase in neutron flux 
above its steady value will be constant with time and reach 


the value: . 
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E x y. Aal super (5.46) 
E uo o Eo) Cols kd: 


where A, is given by equation (5.26) 


k rt 
A, = BE rn o f | J. y, E ‚3 T p EA d 
(U a (Y 1) F (^, E)cesp, E dg dh 


2 2 
Db. ITEN- Ho to 


The zeroth approximation for asymptotic equili- 


brium state is: 


DOR SETS (5.47) 
E. -TIR =Z TE + y? 
O O 

ET 00 ) = 0 (5.48) 


It means that after the transient the neutron flux will 
return to the initial steady value as the stable equilibrium 
state. It is noted that in the absence of delayed neutron 


(8 = 0) and feedback (a = 0), equation (5.46) becomes: 


een en) = AJ, (Mr) cos F (5.49) 
for K = me t y*. 
O O 

This result is exactly the same as that one obtained in 


previous work [6]. 


2. First Approximation 


2 2 
a. INE Ho eya 


The neutron flux will approach the stable equili- 


brium value: 
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Plr, s, =) = D qu 2 _.- = | (5.50) 


as time goes to infinity. < 
2 2 
< 
E. Ho Le 
The equilibrium state of neutron flux is the 


initial steady value 9 (1,2) 
P(r, &,0 ) = $( r, 2) (5.51) 


Therefore, the same physical interpretation can be made 

for a finite cylinder-reactor core as for an infinite 
slab-reactor core. It seems that the asymptotic states for 
a nonlinear nuclear reactor do not depend upon geometric 
configurations. This is true since the transformation is 
independent of coordinate systems and the transformed 


equation is maintained in the same form. 


AL 





VI. CONCLUSIONS AND RFCOMMENDATIONS 


The objective of this thesis has been to provide an 
approach for directly solving the nonlinear space-time 
reactor kinetics equations. The theory is developed for an 
one-velocity, bare, homogeneous, one-group delayed neutron 
and Newtonian cooling system where the temperature rises 
instantaneously with the neutron flux and the temperature 
coefficient is negative (negative prompt feedback model). 

The fundamental tool in this analysis is the Wilhelm's 
nonlinear transformation [8] that transforms a nonlinear 
equation to a weakly nonlinear one in which the nonlinear 
term represents a small perturbation. Then, the reduction 
of the later equation to a linear, parabolic type is estab- 
lished by any method of approximation. The successive 
approximations technique has been used in this thesis since 
it yields a general form of the ith approximation for the 
solution and it is easy to simulate the problem by digital 
sine if need be. 

Furthermore, this transformation is independent of the 
system of coordinates. In other words, the transformed 
equation is maintained in a general form whatever coordinates 
are chosen. Therefore, the same conclusions of stability 
can be drawn from the resulting solutions for reactors of 
different geometric configurations. The asymptotic stable 


equilibrium states are governed by the relationship of the 
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generalized buckling K to u^ (on u^ +77), the lowest eigen- 
value of the associated linear Helmholtz's equation. The 
value, $c (x) + =, is o n state of neutron flux when 
K > u^ (or u^ LTE When K < me the neutron flux returns 
to its initial steady value $o (X) after some transient. 
Then, these results are compared with those derived from 
Theorems I, II and III of Kastenberg and Chambré [1] in the 
case of no-delayed neutron. As mentioned in part IV-E, the 
only inconsistency of the present work with Kastenberg and 
Chambré's is the stable equilibrium state for K = uo. But 
if the comparison is made with linear reactor analysis (no 
delayed neutron and no feedback), the zeroth approximation 
(Cf. see part IV-D-1) will give an asymptotic stable value 
which matches very well with the previous result [6]. So, 
for future work, it is recommended that a verification of 
the first approximation at K = u^ by another approach, say 
digital computers, be made to insure its consistency with 
the theorem III of Kastenberg and with the physical point 
of view. | 

It is noted that the successive approximations CH 
can be applied directly to equation (2.40) to obtain the 
solution provided that the coefficients a and 8* should be 
very small compared to other terms. On the other hand, the 
Wilhelm's transformation reduces equation (2.40) to equation 
(3.10) where the terms n» and No become smaller, even if a 
and B* are not appreciably small. 


Finally, it is assumed from the heginning of this work, 


that the feedback coefficient a is negative. May the 


Ye 





transformation work with a positive feedback system (a > 0 


or a < 0)? Moreover, if Y > 0, u(x,t) should be non-negative 


K 


Kr Otherwise, the func- 


and less than the value, - 
-1) 


tion v$(x,1) will be negative, i.e., the range of u(x,t) is 


a (e 


limited in this case 


o EEN ee x [0 
» Wet. 1) d 


and the singularity may occur when u(x,t) is equal to 


K 


RT at any point (x,T). Therefore, it is not con- 


-1) 
venient to apply this transformation to a problem with 


ale 


positive feedback. 
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APPENDIX A 


—Ó 

UNIFORM CONVERGENCE OF THE SERIES PX e cosu, x cosu,£ 
n=0 

Consider the above series where Un = (2n + 1) 3 . We 


will prove that this series converges uniformly for 


Ea. 45 X, E < iJin 0 < t «T$ e 


O 


L + 2 E «x,t F< a 
-Po T —My T = Mm To E e 
) Lor 
Q) [e "esp veo ELE e es | Tog té% 
2 
Hato 
Expanding e into Taylor's series yields: 
L 
Ha To 4, 2 
e = ENDE RTT)... 
ET 
HaT 
2 
hence, C s > E 
v 
—Kn Te E 
and HL = — —w 
An "Uo Ti To (2n+') 
Therefore, inequality (1) can be written as: 
a) BRT | + ATL | Agr La 
2. e. Cosa xcoc «CIL Nu r 
Pp Mn $ b uL "E (2n 41)” To £t Lo 
It is known that the series X W wets convergent, then 
n NEN = 
inequality (2) implies that the series Y e s cosu,X cosu, 6 
n=0 
should converge uniformly for -a < x, € < a and S. EMO 
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od Er 
Let S(x,8,r) = €. cosm,xcosk, E 


n= 0 


Since each term of the series above is continuous, the func- 
tion S(x,&,T) is also continuous and bounded [11] for 


EN 1 < a and! vac < =. 
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APPENDIX B 


l. Calculation of the Coefficient A, (x) in Equation (4.100) 


Assume K ET. - (sa 


From equation (4.101); 


a a 
A = | A, (x, €) AE + Aa (x, 8) 45 
- 4 sa 
where 
ren 2. 228,2 a tes am ME (1) 
(+A t5) es I+ Z B cospe č 
and 
A (Gt) = Allta) =e Cos M X Cas "ME (1-- —£ ces s E ] (2) 
4 A a ^c. K 0 
a 
Let us evaluate integral /f A. (x,E) dé 
-a 
A. 


= Cos ~ Sin” Ex 
E a Cos 40 Ho H d 
x me 2A Ll Coste, § > 


[ Aso.) de 
Ee 


a T x gt Mo aen) — (3) 


+ E 6, cos pto E 

Let D NE == (4) 
TE > 
dy = -Hosinpz de 


Equation (3) becomes: 
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a 


l pp 
[5e 04 == Resah | E K J 


The integral in the right hand side of equation (5) is 


separated into 2 different integrals as follows: 


= [VEE dy- p| EE dy “ 


But the first integral in the right hand side of equation 


(6) can be calculated from C.R.C. tables, 


Kar: dy = E (7) 


The second integral can he evaluated after changing of 


Variable: : 
Let 2 52 ww Y + p 

Ec eu 
hence, 


= ^ 
E Vi- G dy E | V-xt-2py«ni- v? m (8) 
o p > 


Since it is obtained from any integral table that: 
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are dE VT Spe LI P + of Ax 
x: umo ZEAL p* 


| er (9) 
Ya pl pr 


dx ate er 
and T z =~ sin un SE (X-p) 
V AAE K l p* Va de P 








i z sin E (10) 
An | Ç "( X e| — =) 
and a ein f P 
AE p=! S 
= l sn” ERL (11) 
Vpr- y+p 


Substitute equation (10) and equation (11) into equation (9) 


to obtain: 
| E +1- pt = E o wor Sa = i^ 
p we i-pt de lVIDQT psu Np sin Cd (12) 


which is integral (8). Thus we have 


i : 
| Vi- y" dy = — W V pi +Vp=1 ein os (13) 
EE nr P 


Since p = =a is very large, hence, 
O 
A 
P 
V == N z ne E aan. s 
and p = bL 4.) Tire a 2p gp> 
c MESA 
ey) — ae 
P I? cp” 


Equation (13) can be written as follows: 
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i 


WT 
f <— dy ag oo T(p-a- 2p T +(p-4 DE i» n d ) 


ger 
(14) 


© 


qe. 


IA 


"BTR 
Now, inserting equation (7) and equation (14) into equation 


(6) yields: 
«Bb 


Zu = Te = E 
Be Ay J = Phas - so = |. = 3K 


id 
which is the value of the integral in the right hand side of 





equation (5) Thus, 


å 

A A Eg S 
lar 5 S 

m 
But Ro = > 
Then 
a 
*- 
| A CAS: = = EX ES CoS MaX (15) 
a 
(x,E)dE where A,(x,8) is given 


Next, evaluate integral f A 
| -à 


by equation (2) 

E Y a 
ij A ax, x) Ax = i 8, coses | (14 becos Eyes put de 
-a AÀ e Fi K | 


A a 
o COS KL pot Cos" Mey Ax +28, cos ds & te] $ é ) 


£g 


Ode 
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dy = yz de 
hence, 


M 


A- a Ex 
p Uu 7. 
f Cos MoE dg E Cos Ho 5 dz E E Cos y 


E z (17) 
ELSE + 


and similarly, 


— 


a It 
3 L- 
J cos tes dg = MT cst dye e hh) 
o 


= ee (18) 
3 TC 


Substituting equations (17) and (18) into equation (16) gives 


N | 
E 
[ AaGizdde = E B,cen xla+8 TB. 


hi 





P P ae ai 
C -on 


Thus, A, (x) is obtained by summing up equations (15) and (19) 


Ae = -5 I, e Bo cos by x + (E + E Estas P, costo x 
a. 
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or 


= {6° & BK Pb, 
Ag = (£44 Gets LE RIR cosp (20) 


Ac E EIG 


Let us consider the order of magnitude of each term in brac- 


kets. To do this, we assume the following data [13] 





A z 5 ß = ©. 00 323 
> = n Zo o 236 DF = 0.0052 em! 
< 

C == IG a fom Cp = 0.07 E cal [9% 

7 -5 © zu 
V z u-?Glxlo' cr [sc Xy —S5xlo K 

—| =) 
A = 0-824 sec e = 200M4Y = 200n3.82n06 col 
fission Lission 


Therefore, we obtain: 


2 
“=(E) = (5) = 00986 


K o = 7? S = I. 0986 


s -l6 
NEP END UTE SO ma a Ae. 
Verde: 
T 2 - 
= [Li-@) xo, | = | AY) = S.S xio Cao e. C 
A l 6 
^c E ——— MS 2329 = 32.2 «l0 | C di nia Sion less ) 
E. 4. Sélxto! x 0-005 t 
— 20 
BoA = 4..7 x10 dU LN : 
3 TK 
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i or uc cue | (21) 


Â x $ -21. 2. 
3 Pcr x 2 cte (22) 


«|n 


T 
A 


mum 
mS = #.6xI0 cm” oec. (23) 


g^ = 0.06236 AÀ«. = 0.00 362 4 x16 ° = UU x10 


Thus, comparing equation (21) to equations (22) and (23) 


yields: 
ye — 
Ep MESSI F 
3 WK AG A 
and Sae E 6> 
Jat ow bo 


and the value of A, (x) is approximately equal to 


X 
~ 24 
Kh e ER becos Ey (24) 


Hence, A.(x) is finite and positive for -a < x < a. 
2. Calculation of the Coefficient Ag (x) in Fquation (4.114) 
2 


Assume K > Ho 


From equation (4.115): 


A 9 (x) 


1! 


[Arc S) dg 
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where A che ie De Q «x Bo 


— 


Ces X Co£ ( S 
O t) PS P" 


o 


a 
= 
BR) - Rea B _ ces x J| cos "c „E)d 
L AKALLK- Kor) T -— 4 à 5 er 


Since the value of the integral in the right hand side of 


expression (25) is given by equation (18), hence, 


= 
X 
Aa (x) = 5 _ f*xBo_ coc x 
a 


(26) 
NKA CK- (odo 


The equation above shows that A, (x) is finite and positive 


tor -a < X < a. 
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